The Functional Integration and 
'"§'■ the Two-Point Correlation Function of 

(n: the One-Dimensional Bose Gas 

9^' in the Harmonic Potential 

oo 



m 



c^ 



X 



N. M. Bogoliubov*, C. Malyshev^ 

V. A. Steklov Institute of Mathematics, 

St. -Petersburg Department, 

Fontanka 27, St. -Petersburg, 191023, RUSSIA 






-)— > 

a 

O ' Abstract 

o : 

A quantum field-theoretical model which describes spatially non-homogeneous one- 
dimensional non-relativistic repulsive Bose gas in an external harmonic potential is 
1^ ' considered. We calculate the two-point thermal correlation function of the Bose gas in 

the framework of the functional integration approach. The calculations are done in the 
Tj- i coordinate representation. A method of successive integration over the "high-energy" 

^ I functional variables first and then over the "low-energy" ones is used. The effective 

\^ ■ action functional for the low-energy variables is calculated in one loop approximation. 

^^ , The functional integral representation for the correlation function is obtained in terms 

of the low-energy variables, and is estimated by means of the stationary phase approx- 
imation. The asymptotics of the correlation function is studied in the limit when the 
temperature is going to zero while the volume occupied by non-homogeneous Bose gas 
'^ I infinitely increases. It is demonstrated that the behaviour of the thermal correlation 

" ■ function in the limit described is power-like, and it is governed by the critical exponent 



O 

O ■ which depends on the spatial and thermal arguments. 
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§1. Introduction 

A recent burst of interest to theory of the Bose gas is caused by experimental realization 
of Bose condensation in the ultra-cold vapors of alcali metals confined in the magneto- 
optical traps [1,2]. In particular, it became possible to study the Bose condensation in 
the systems, which are effectively two-dimensional or quasi one-dimensional [3,4]. Here a 
partial localization along one or two directions in three-dimensional system is achieved by 
making the level spacing of the trapping potential in the corresponding directions larger than 
the energies of individual atoms. The field models which describe the Bose particles with 
the delta-like inter-particle coupling confined by an external harmonic potential provide a 
good approximation for a theoretical approach to the experimental situation [5]. The theory 
of non-ideal Bose gas attracts traditionally not only physicists but also mathematicians 
[6-8]. For a translationally invariant homogeneous Bose gas, the field model in question 
corresponds to a quantum nonlinear Schrodinger equation, which admits an exact solution 
in the one-dimensional case [9, 10]. This fact allows to obtain closed expressions for the 
correlation functions [11,12]. 

In real physical systems the interest to a transition from three-dimensional to one- 
dimensional behaviour is caused by the fact that an effective density of atoms in the one- 
dimensional Bose gas can be either high or low depending on the parameters of the sys- 
tem [13]. Here a low density implies a strong coupling between the particles [14,15], while 
a high density corresponds to a weak interaction. 

The present paper continues the series of the papers [16-20] devoted to investigation of 
the correlation functions of the Bose gas with weak repulsive inter-particle coupling in the 
presence of an external harmonic potential. Since there is no exact solutions in the case of an 
external potential, in the present paper, as well as in Refs. [16-20], the method of functional 
integration [21-24] is accepted for investigation of the correlation functions. In the present 
paper we develop the results of the paper [20], where, in a distinction with [16-19], an explicit 
dependence of the correlation functions on the imaginary time for non-zero temperature was 
taken into account. It will be demonstrated that the presence of the external potential 
(of the trap) results in a change of the asymptotical behaviour of the two-point correlation 
function in comparison to a translationally invariant case. It will be demonstrated that the 
observed change happens in the range of temperatures, which are comparable with inverse 
of the characteristic length of the trap provided this length tends to infinity. 

The paper is organized as follows. Section §1 has an introductory character. A descrip- 
tion of the one-dimensional model of non-relativistic Bose field in question, as well as a 
summary of the method of the functional integration, are given in Section §2. An approach 
to approximate investigation of the functional integrals is also presented in Section §2. This 
approach is based on a successive integration first over the high-energy over-condensate ex- 
citations and then over the variables, which correspond to the low-energy quasi-particles. 
Besides, in this Section we give a derivation of one loop effective action for the low excited 
quasi-condensate fields, and the corresponding energy spectrum of the low lying excitations 
is obtained. The method of stationary phase is used in §3 for an estimation of the functional 
integral, which expresses the two-point thermal correlation function of the non-homogeneous 
Bose gas. The method of asymptotical estimation of the correlators, which is used in the 
present paper, was proposed in [25], where the asymptotical behaviour of two-point Green 
functions of the homogeneous Bose gas was investigated for the spatial dimensionalities one, 
two, and three. In the present paper it is demonstrated that the stationary phase method [25] 
admits a generalization for the spatially non-homogeneous Bose gas in the external potential 



also. The asymptotics of two-point correlation functions of non-homogeneous Bose gas are 
obtained in Section §4 both at nonzero and zero temperatures. A short discussion of the 
results of the present paper is given in Section §5. 

§2. The effective action and the Thomas-Fermi approxi- 
mation 

1. The partition function 

In the present paper we shall consider one-dimensional Bose gas described by the Hamiltonian 
H defined on the real axis M 9 x: 

/ I '^''" (^) '^'^ (^) "I — "^^ (-'') ^"^ ('*') "^ (•^) "^ (-'') I '^•'^^ 

where ^"^(x) and t/'(x) are operator-valued fields, which describe creation and annihilation 
of the quasi-particles over the Fock vacuum |0). The fields ^'^(x) and ^(x) are subjected to 
the commutation relation 

^{x)^''{x') — ■^^(a:')-^(x) = 5{x — x') 

(i^'^ix) and ^(x) are mutually commutative), and H is the "single-particle" Hamiltonian. 
The following notations are used in equation (1): m is the mass of the Bose particles, /x is 
the chemical potential, g is the coupling constant which corresponds to the weak repulsion 
(i. e., g > 0), and the external confining potential is chosen in the form of the harmonic 
potential V{x) = fn'^x\ 

Let us begin with the investigation of the partition function Z. It can be represented in 
the form of the functional integral [21-24]: 

e'^l^'^lpV'^V', (2) 

where S[ip,i/j] is the action functional of the system in question: 

d 



S[iIj,iIj]= dr dxiilj{x,T)(— H]iP{x,t) 



V'(x,r)V^(x,r)V^(x,r)V^(x,r)j. (3) 



The domain of the functional integration in (2) is given by the space of complex-valued 
functions ijj{x,t), tp^x^r) depending on two arguments: x G R and r G [0,/3]. With regard 
to the first argument x, the functions ^(x, r), ^(x, r) belong to the space of quadratically 
integrable functions -L2(R), while they are finite and periodic with the period 13 = {ksT)'^ 
with regard to the imaginary time r (ks is the Boltzmann constant, and T is an absolute 
temperature). The variables ^, ^ are the independent variables of the functional integration 
[21], and VipVip is the functional integration measure. 



At sufficiently low temperatures we can expect to assume that each of the variables 

■?/'(a;,r), ■?/'(a;, r) is given by two constituents. One of them, ■?/'o(a:, r), ipo{x,T), corresponds to 
quasi-condensate, while another one — to the high-energy thermal (i. e., over-condensate) 
excitations '?/'e(a:, r), ■?/'e(x, r): 

^(a;,r) =^o(x,r) +^e(a:,r), ^{x,t) = ^o{x,r) + ipeix^r). (4) 

It should be stressed that here and below just the quasi-condensate is assumed, since a true 
Bose condensate does not exist in one-dimensional system [6]. In the exactly solvable case, 
the existence of the quasi-condensate implies that a non-trivial vacuum state (i.e., a ground 
state) exists. The quasi-condensate variables ^o(a;, t), z/'o(x, r) can also be represented as 
the sums of the constituents: 

iJo{x,T) =^po{x) +^{x,t), i)o{x,T) =2po{x) +^{x,t), (5) 

where the field ipoix) describes the ground state of the model at zero temperature, while the 
field ^{x,t) describes the low lying excited particles. Let us require the variables in (4) to 
be orthogonal in the following sense: 

ll)o{x,T)ll)e{x,T)dx = I tl>o{x,T)^lJeix,T)dx = 0. 

As a result, the integration measure VipVip will be replaced by Vipo'I^'4'o'^'4'e'^'4'e- 

To investigate the functional integral (2), we shall perform a successive integration over 
the fields ^,^. First, we shall integrate over the high-energy constituents, and then - over 
the low-energy ones (see. (4)) [21,24]. At a second step, it is preferable to pass to new 
functional variables, which describe an observable "low-energy" physics [21,24] in a more 
adequate way. After the substitution of the expansion (4) into the action (3) we shall take 
into account in S only the terms up to quadratic in ipe, ipe- This means that we are making an 
approximation in which the over-condensate quasi-particles do not couple with each other. 
In this case, it is possible to integrate over the thermal fluctuations ^e(a^, t), ipe{x,T) in a 
closed form and thus to arrive to an effective action functional 5'eff [t/'o,^o] depending only 
on the quasi-condensate variables ipoyipo'- 



where the tilde in S implies that "self-coupling" of the fields ilJe{x,T), ipe{x,T) is excluded. 
With respect to equation (6), the partition function Z of the model takes an approximate 
form: 

Z^ /e^^«l^°'^°'P^«PV'o. (7) 



Let us consider the derivation of the effective action Sefi[ipo,'4'o] (6) in more details. The 
splitting (4) allows to derive Ses[ipo,4'o] in the framework of the field-theoretical approach 
of the loop expansion [26,27]. We substitute (4) into the initial action S'[^,^] (3) and then 
pass from S to the action S, which is given by three terms: 

•J = 'Jcond + '-'free + '-'int- (oj 



In (8), S'cond is the action functional of the condensate quasi-particles, which corresponds to 
a tree approximation [26,27]: 



g - - "1 

- -1po{x,T)^o{x,T)iJo{x,T)iJo{x,T)j. (9) 

At the chosen approximation, the action for the over-condensate excitations 5'free takes the 
form: 

^frce[^e, i^e] = \j dr j dx{^Pe, ^e)^"^ ("^A . (10) 



Eventually, ^int is given by the part of the total action functional, which describes the 
coupling of the quasi-condensate to the over-condensate excitations: 



+ llJe{x,T)[K^ - giljo'ipo\i>o{x,T)\. (11) 

In formulas (9)-(ll), we have defined the differential operators K± = ±d/dr — Ti (here the 
Hamiltonian Ti is defined in (1)) and the matrix-differential operator G^^: 

G-' ^ G,' - % (12) 

where 

In this approach it is appropriate to apply the stationary phase method to the functional 
integral (6). To this end, let us choose ^o, ipo as the stationarity points of the functional S'cond 
(9), which are defined by the extremum condition (5 ( S'cond [^o,^o]) = 0. The corresponding 
equations are taking the form of the Gross-Pitaevskii-type equations [28]: 



d h^ & 



Or 2m dx^ 

d h^ &■ 



+ /i - v{x)j^o - g{^o^o)^o = 0, 



. _ _ (13) 

Or ' 2mdx^^^~^^^'^)^°~^^^°^°'^^°^^' 



The contribution of the action S'int (H) drops out from (8) since ipoy'ipo are chosen to be 
solutions of equations (13). Therefore the dynamics of ^e,^e is described, in the leading 
approximation, by the action S'free (10). The latter depends on ipo,'ipo non-trivially through 
the matrix of the self-energy parts S, which enters into G~^ (12). The Thomas-Fermi ap- 
proximation is essentially used in the present paper in order to determine the stationarity 
points ipoti'o- This approximation consists in neglection of the kinetic term 2^^ in equa- 
tions (13) [5,28]. The Thomas-Fermi approximation is valid for the systems containing a 
sufficiently large number of particles, and it is widely used in the theoretical approaches 



to description of the Bose condensation in the magneto-optical traps [5,28]. The follow- 
ing condensate solution can be obtained provided only r-independent solutions of (13) are 
allowed: 

'ipo'ipo = Ptf{x; iJ.) = -{fi- V{x)) 6 (/i - V{x)) , (14) 
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where 6 is the Heavyside function. Now the integration in (6) with respect to ipe, 4'e is 
Gaussian. This leads [21] to the one loop effective action, which takes the following form in 
terms of the variables ipo, ^o'- 

^eff[^o,^o] = 5cond[^o,^o] - ^ In Det(G-i). (15) 

Here G~^ is the matrix operator (12), and ipo, ipo have a sense of the new variables, and 
their dynamics is governed by the action (15). 

In order to assign a meaning to the final expression for the effective action (15), it is 
necessary to regularize the determinant Det(G'"^). In our case, the operator G~^ is already 
written as a 2 x 2-matrix Dyson equation (12), where the entries of S[^o,^o] play the role 
of the normal (Sn = E22) and anomalous (S12, S21) self-energy parts. The Dyson equation 
(12) defines the matrix G, where the entries have the meaning of the Green functions of the 
fields ipc'ipe- The matrix G arises as a formal inverse of the operator G~^: 

g=(g,'-y:Y\ (16) 



The matrix operator G~^ (12) can formally be diagonalized by means of the famous N.N.Bo- 
goliubov's (m, w)-transform [6]. The corresponding equations, which describe the unknown 
coefficient-functions u, v, result in a compatibility requirement, which determines, in turn, 
the quasi-classical spectrum of the elementary excitations [28]. 

With regard to our purposes, it is appropriate to represent G~^ as follows: 

G-' = Go' -^ = g-'-{^- 2gpTF{x;fM)T), (17) 

where / is the unit matrix of the size 2x2, and the matrix Q~^ is defined as 

K+ - 2gpTF{x- fi) \ _ flC+ 



K_-2gpTFix;fi)J ~ \0 /C_ 

Here ptf{x; fj,) is the solution (14), and equation (17) implies that we simply added and 
subtracted 2gpTF{x; p) on the principle diagonal of the matrix operator G~^. A formal 
inverse of the operator Q^^ can be found from the following equation, which defines the 
Green functions Q±: 

Using the relation In Det = Tr In, one gets: 
-ilnDetG-i = -iTr In {? - g{^ ~ 2gpTF{x- p)l)^ - ^InDet (^^ ^). (19) 



The first term in Right Hand Side of (19) is free from divergencies. Let us consider the 
determinant of the matrix-differential operator in Right Hand Side of (19). The operators 
K,± can be written in the form: 

Let us denote the eigenvalues of the operators /C± as ±.iuB — A„, where ujb are the bosonic 
Matsubara frequencies, and A„ are the energy levels (which are labeled by the multi-index n) 
of the operator — ^^ — \y{x) — /u|. The regular part of the logarithm of the determinant 
in equation (19) has the sense of the free energy F„c of the ideal gas of the over-condensate 
excitations: 

where the regularized values of the determinants of the operators /C± can be obtained, for 
instance, by means of zeta-regularization approach [27]. Then, in the leading order in g, one 

gets: 

- -In Det G-^ 
2 

P 
^ -/3F„c(/i) + g dr dx{g+{x, r; x, r) + G-ix, r; x, r)) {ipoipo - Ptf{x; /i)) 



= -(3Fnc{n) +2g dr dxpnc{x)i)oi)o- (21) 



Here Fnc is the free energy of the non-ideal gas of the over-condensate quasi-particles, and 
the last term in (21) describes a coupling of the over-condensate quasi-particles with the 
condensate. The density of the over-condensate quasi-particles is Pnc{x) = —Q±{x,t;x,t), 
and it depends only on the spatial coordinate x. At very low temperatures and sufficiently 
far from the boundary of the domain occupied by the condensate, the quantity Pnc{x) can ap- 
proximately be replaced by Pnc(O), since Q±{x, r; x, r) is almost constant over a considerable 
part of the condensate [29]. Eventually, one obtains 

'S'eflf[V'o,V'o] = -pFncip) + dr dx^iJo{x,T){^—-—^ + A-V{x)jtpo{x,T) 



-^i'o{x, T)i!o{x, T)^poix, t)^Po{x, r) J, (22) 

where A = fi — 2gpnc{0) is the renormalized chemical potential. We shall consider Ses (22) as 
one loop effective action, where the thermal corrections over the "classical" background are 
taken into account. The "classical" background corresponds to the solution (14). It should 
be noticed that the described derivation of the effective action does not depend on spatial 
dimensionality, and therefore it is valid for two and three dimensions also. 

In the effective action (22) it is appropriate to pass to the new variables, namely the 
density p{x,t) and the phase v^(x, r) of the field ^o{x,t) [21]: 



M^, r) = V;K^e^^^"•"^ M^, r) = ^/^(^)e~'^^^^^\ (23) 



We shall consider p(x, r) and ip{x^ r) as two new independent real-valued variables of the 
functional integration. Now the integration measure Vipo'Dipo is replaced by the measure 
VpVif. In these new variables the effective action (22) takes the form: 



+ I'dr I dx{^{^dl^-p{d^Vr) + {A-V)p-^-p'}. (24) 

Here and below we denote the partial derivatives of the first order over r and x as dr and 
dx, respectively, whereas the partial derivatives of the second order — as d^ and d^. Notice 
that equations (22) and (24) remain correct at K = 0, so that the renormalized chemical 
potential is still given by the equation A = p — 2gpnc{0), where Pnc(O) is the "bare" [21] 
density of the condensate. 

2. The excitations spectrum 

Let us consider the problem of determination of the spectrum of the low-energy quasi- 
particles. We shall apply the stationary phase approximation to the integral (7), where 
the effective action is given by (22). The corresponding stationarity point is determined 
from the extremum condition S(^Sefi[p,ip\j = 0, which is equivalent to the couple of the 
Gross-Pitaevskii equations: 

^dr^ + ^ i^dl^p - {dxvf) + A - V{x) -gp = 0, 

-idrp H dx {pdx(p) = 0. 

m 

Let po and ipo to denote some solutions of the couple of equations (25). Substituting po,ipo 
into the effective action (24), one obtains 

Ses[pO,fo] = -pFnc{p) + o / ^^ / ^^'^O' ^'^^^ 



Here Fndp) is the free energy of the non-ideal gas of the over-condensate quasi-particles. 
The total free energy of the system is F{p) = —^Ses[po, V^o] [18]. 

Let us use the Thomas-Fermi approximation, which is valid at sufficiently low tempera- 
tures, and drop out the kinetic term (c^^\/p)/a/p in the first equation in (25). Solution with 
drP = = drf appears, provided the velocity field v = m'^dxf in (25) is taken equal to 
zero. In this case, equations (25) lead to the density of the condensate described by the 
solution (14), where the chemical potential p is replaced by A: 

^.W.^>.W = ^(i-J)e(i-J). (27) 

Explicit form of the external potential V{x) = yfi^x^ is taken into account in expression 
(27). The form of the solution (27) means that the quasi-condensate occupies the domain 
|x| < Re at zero temperature. The length Re defines the boundary of this domain, R^ = 
:^2 (in three dimensional space, this would correspond to a spherical distribution of the 



condensate). In the homogeneous case given by the hmit 1/Rc -^ 0, the Thomas-Fermi 
solution Ptf{x) is transformed into the density Ptf(O) = A/g, which coincides with the 
density of the homogeneous Bose gas [21]. 

Following the initial splitting (5), we suppose that the thermal fluctuations in vicinity 
of the stationarity point (27) are small, and therefore an analogous splitting can be written 
for the condensate density also: 

poix, t) = ptf{x) + tiq{x, t). (28) 

The Gross-Pitaevskii equations (25) linearized in a vicinity of the equilibrium solution ptf = 

Ptf{x), ^ = const, then takes the form: 

idr'^Q - gno + S^ttq = 0, 

^mpTF (29) 

idriiQ dx{pTFdx(po) = 0. 

m 

Eliminating (po and dropping out the terms proportional to K^, we go over from (29) to the 
Stringari thermal equation [30]: 

^a% + 9.((l - ^)5.vro) = 0, (30) 

where the parameter v has a meaning of the sound velocity in the center of the trap: 

v' ^ P-I^^ = A. (31) 

m m 

The substitution ttq = e'^'^'^u{x) transforms (30) into the Legendre equation: 

''\u{x) + ^((l-^)^u{x))=0. (32) 



h?v'^ dx\\ R^J dx 

Since the Thomas-Fermi solution (27) is different from zero only at |a;| < Re, we shall 
consider equation (32) at x G [—Re, Re] C M, as well. After an analytical continuation 
uj — > iE, equation (32) possesses the polynomial solutions, which are given by the Legendre 
polynomials Pn{x/Rc), if and only if 

^)^^'-sIf^= ="("+!). "SO- (33) 

In other words, equation (32) leads to the spectrum of the low lying excitations: En = 



MlJ ^^'^^ ' , n > [31]. Notice that the corresponding equation for the homogeneous Bose 
gas is obtained after a formal limit l/Rc — > in (32) at finite x. Provided the latter is 
still considered for the segment [—Re, Re] 3 x with a periodic boundary condition for x, 
we arrive at the discrete spectrum of the following form: Ek = hvk, where k is the wave 
number, k = {7i/Rc)n, n E Z. 



§3. The two-point thermal correlation function 

Let us go over to the main problem of the present paper — to the calculation of the two- 
point thermal correlation function of spatially non-homogeneous Bose gas described by the 
Hamiltonian (1): 

T{xi,Ti;X2,T2) = {Tr^\xi,Ti)i}{x2,T2)) , (34) 

where Tr is a "r-chronological" ordering with respect of the imaginary time r, and the 
angular brackets (, ) correspond to averaging with respect of the Gibbs distribution [6]. We 
can express the correlator r(xi, ri; X2, T2) as a ratio of two functional integrals [21-24]: 

r(a;i, n; x^, r^) J^^m^^^ ' ^^S) 

where the action S'[^,^] is defined in (3). 

We are interested in the behaviour of the correlators at the distances considerably smaller 
in comparison with the size of the whole domain occupied by the condensate. The main 
contribution to the behaviour of the correlation functions is due to the low lying excitations 
at sufficiently low temperatures [21]. To calculate r(a:i, ri; 0:2, T2) (35), we use the method of 
successive functional integration first over the high-energy excitations ^e, ^e, and then over 
the low-energy excitations V'o, V'o [21,24,25]. We find that in the leading approximation the 
correlator we are interested in looks as follows [21,25]: 

^, . Je'^^«^^'"^''^Mxi,n)Mx2,T2)v^oV^iJo f^^. 

^^-----^^-^) Je^e.[^^^^^lP^.PV^. ' ^^'^ 

where 5'eff[^o, ^o] is the effective action (22). We can rewrite (36) in terms of the density- 
phase variables (23), and then represent the integrand in the nominator in the form of a 
single exponential: 



r(xi,ri;x2,r2) 

/exp (S'cflf[p, v?]-iv^(a;i, Ti)+iip{x2, T2) + \ lnp(a;i, ri) + | lnp(x2, T2)]VpVyD 

/exp {Scs[p,'^])T>p'Dip 



(37) 



Here Ses [p, <^] is the effective action (24) . 

Since the fluctuations of the density are suppressed at sufficiently low temperatures [28], 
one can replace lnp(xi,ri), lnp(a;2,T2) in (37) by InpTpixi), lnpTi?(x2), where the density 
Ptf is defined by (27). In accordance with the variational principle suggested in [25], we 
shall estimate the functional integrals in (37) by the stationary phase method, and we shall 
consider only a leading approximation. For the correlation function T{xi,Ti;x2,T2), we 
obtain the following estimation: 



r(a;i, Ti; X2, T2)'=:i\/ Ptf{xi)Ptf{x2) 

X exp(-S'cflf [po, v^o]+'S'cff [Pi, ^i\-ivi{xi, Ti)+iipi{x2, ts)), (38) 

where the variables Po,<fo are defined by the extremum condition 6{Ses[p,<^]) = 0, and 
therefore they just satisfy the stationary Gross-Pitaevskii equations (25). For Ses[po, 'Po] we 
use the expression (26). 
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The fields pi,^i are defined by the extremum condition: 

6(Scs[p, (p] - i(p{xi, Ti) + i(p{x2, T2)] = 0. (39) 

This variational equation leads to another couple of equations of the Gross-Pitaevskii type. 
One of these equations turns out to be a non-homogeneous equation with the 5-like source, 
while another one — a homogeneous equation. In fact, the homogeneous equation appears 
due to a requirement of vanishing of the coefficient at the variation 6p{x, r), 

^^rV + ^ {^dl^ - (d.^y) + A - Vix) -gp = 0. (40) 

In its turn, the non-homogeneous equation is defined by vanishing of the coefficient at the 
variation S(p{x,T): 

—idrP H dx {pdx^) = iS{x — xi)6{t — ri) — i6{x — X2)S{t — T2). (41) 

Substituting the solutions pijV^i, which respect (40), (41), into the effective action (24), one 
gets 

Sesipi, ^1] = -f3Fncip) -^ + ^[drf dxpl (42) 



Further, it can consistently be assumed that the solution pi(x, r) can be represented 
as a sum of Ptf{x) and of a weakly fluctuating part provided the boundary Re is far from 
beginning of coordinates: pi(a;, r) = pTF{.x)+ni{x,T) (for a comparison, see (28)). Therefore, 
the terms a/^JTc^xa/^ and dxTiid^fi are small and can be omitted. Taking into account a 
linearization near the Thomas-Fermi solution, one can finally re-write equations (40) and 
(41) as a couple of the following equations: 

idr(fi - fi-TTi - -— {dxfif = 0, (43.1) 

2m 

—idrT^i H dx (prrdx'^i) = iS(x — xi)6(t — ti) — iS(x — X2)S(t — T2). (43.2) 

m 

Differentiating (43.1) over r, substituting the result into (43.2), and dropping out the 
terms of higher orders in g and /i^, one obtains the following equation: 

1 h^ 

-d^ipi H dx {pTF{x)dx(pi) = iS{x - xi)6{t - n) - i6{x - X2)S{t - T2). (44) 

g m 

It is convenient to rewrite it as follows: 

1 TflQ ( '\ 

dlipi+dx {'pTF{.x)dx'{)i)=i^::r-{ 5{x - xi)5{t-ti)-5{x-X2)5{t-T2) k (45) 



where v means the sound velocity in the center of the trap (31), and ptf is defined by (27). 
More specifically, solutions of equations (44), (45) depend on the coordinates of the 5-like 
sources in Right Hand Side, i.e., on Xi, ri, X2, T2: ipi{x,T) = ipi{x,T;xi,Ti,X2,T2). Now, 
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with the help of equations (26), (42) and (43), one can calculate the contribution, which is 
a part in the exponent in (38): 

- <S'cflf[po, V'o] + Scs[pi, Vi] 

13 13 

~ 2 1^ c?x(p? - pI) = - dr dx{pi - po)(pi + Po) 



13 13 ^ 

c:^g dr dxTTipo = dr dx\idripi - ^(^xV^i) JPo 







/3 (3 

-.2 





= ^{vi{xi,ri)-'^i{x2,T2)). (46) 

Substituting (46) into (38), one obtains the following approximate formula for the correlator: 

r(xi,ri;x2,r2) ~ \/ptf{xi)ptf{x2) exp \^--[ipi{xi,Ti) - '^i{x2,T2))j. (47) 

It is natural to represent the solutions of equations (44), (45) in terms of the solution 
G{x, t; x', t') of the equation 

^2^lG{x, r; x\ r')+4 ((l-^)9.G(x, r; x', r'))=^^5{x-x')5{r-r'). (48) 

Bearing in mind the homogeneous equation (30), we shall call (48) as non-homogeneous 
Stringari equation. As a result, the representation (47) can be re-written as follows: 

T{xi,Ti]X2,T2) 

- \/ptf{xi)ptf{x2) exp(-- {G{xi, Ti;x2, T2) + G{x2,T2;xi,Ti)) 

+ -G{xi, n; xi, Ti) + -G{x2, T2; X2, Ta) j . (49) 

As it is clear after [8], the function G{xi, Ti, X2, T2) has a meaning of the correlation function 
of the phases: 

G{Xi,Ti;X2,T2) = -{(^{xi,Ti)ip{x2,T2)), (50) 

where the angle brackets in Right Hand Side should be understood as an averaging with 
respect to the weighted measure VpVip exp (S'egfp, v^]). Substituting (50) into (49), we 
obtain the known approximate formula for the correlator (34) [8,21]: 

{Tr^\xi,Ti)ij{x2,T2)) ~ \/ Ptf{xi)ptf{x2) exp [--{{'^{Xi.Ti) - ^(X2, Ta))^) j . (51) 

Notice that the terms in the exponent in (49) have different meanings. The terms 
G{xi,Ti\ X2, T2) and G{x2, T2; Xi,Ti) depend on the differences of the arguments of two-point 
correlation function and therefore are responsible for the behaviour of the correlator at large 
distances. The terms G{xi,ti;xi,ti), G{x2,T2;x2,T2) each depend only on a single set of the 
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coordinates and thus contribute to the amplitudes only. The Green function G(x, r; x', r') 
depends, in fact, on the difference of r and r' due to the invariance under the shifts of r 
(see. (48)). Therefore, only a dependence on the spatial coordinates remains provided the 
corresponding thermal arguments coincide. Then, it is possible to represent the correlation 
function as follows: 

r(a;i,ri;a;2,r2) ~ \/'p{xi)'p{x2) exp y--{G{xi,Ti\X2,T2) + ^(xs, Ts; Xi, ri))j, (52) 

where p(a;i), ]){x2) are the renormalized densities [18]. The solution G{xi,ti^X2,T2) of equa- 
tion (48) is defined up to a purely imaginary additive constant, which has a meaning of a 
global phase. As it is seen from (51), this constant does not influence the fluctuations. 

§4. The asymptotics of the correlation function 

Therefore, the problem concerning the study of the asymptotical behaviour of the two-point 
thermal correlation function r(a;i, ri; 0:2, T2) (34), which is given by the representation (37), 
is reduced to solution of the non-homogeneous Stringari equation (48). The corresponding 
answer (or its asymptotics) should be subsequently substituted into (52). In the present 
section, we shall obtain explicitly solutions of equation (48), and we shall consider the 
corresponding representations for the asymptotics of r(a;i, ri; a;2, r2). Let us begin with the 
limiting case of equation (48), which corresponds to a homogeneous Bose gas. 

1. The homogeneous Bose gas 

In this subsection we shall consider the asymptotical behaviour of the correlation function of 
the homogeneous Bose gas. As it was mentioned above, the homogeneous case corresponds 
to V{x) = 0, and the related equation appears from (48) at 1/ Re — > 0: 

^ -,dlG{x, r; x', r') + dlG{x, r; x', r') = ^6{x - x')6{t - t'). (53) 



We consider (53) for the domain [—Rc^Rc] x [0,/9] 3 {x,t) with the periodic boundary 
conditions for each variable (in other words, we consider (53) on the torus S^ x S^ 3 (x, r)). 
The (5-functions in Right Hand Side of (53) are treated as the corresponding periodic 6- 
functions. This allows us to represent the solution of this equation as the formal double 
Fourier series: 

„ 1 Aui(T—T')+ik(x—x') 



G(x,r;x',r')=(^)(2/3i?, 



fl2^2)y ^ ''^) A^ ^2/(;^2^2) + p 
ui,k "- 

where u = (27r//3)/, / G Z. The notation for the energy E/. = hvk, where k = {7i/Rc)n, n E Z, 
is used in the representation (54). Besides, the representation (54) requires a regularization, 
which consists in neglection of the term given hj 00 = k = 0. 

Using (54), one can deduce two important asymptotical representations for the Green 
function. In the limit of zero temperature and of infinite size of the domain occupied by the 
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Bose gas, one can go over to the asymptotics of r(a;i, ri; X2, r2). When a strong inequahty 
P~^ = ksT ^ hv/Rc is valid, we obtain: 



G{x, r; x', r') 
9 



2Trhv 



ln<^2 



sinh 



vr 



hpv 



\x — x'\ + ihv{T — t')) 



9 






l\2 



ApR^ h'^v^ 



+ C, 



(55) 



where \x — x'\ < 2Rc, \t — t'\ < P, and C is some constant, which is not written expUcitly. 
When an opposite inequahty P~^ = ksT <^ hv/Rc is vahd, we obtain: 



G{x,t; x', r' 

9 

2-71 hv 



ln<^2 



sinh 



tTT 



2Rr 



\x — x'\ + ihv{T — t')) 



APRc 



\t — T \ + C 



(56) 



where \x — x'\ < 2Rc, \t — t'\ < j3, and C is another constant. 

Let us substitute the estimate (55) into the representation (52) and take simultaneously 
the limit phv/Rc -^ (the size is growing faster than inverse temperature). Then, we obtain 
the following expression for the correlator in question: 



r(xi,ri;x2,r2) ~ ^/p{xi)p{x2) 



sinh 



TT 



hpv 



\Xl - X2\ +ihv{Ti - T2)) 



-g/2'Khv 



(57) 



Further, applying the relation (56) and taking the limit Rc/{Phv) -^ (the inverse temper- 
ature grows faster than the size), we obtain for r(a;i,ri; 0:2, r2): 



T{xi,Ti;X2,T2) ~ ^/p{Xi)p{x2] 



sinh 



2Rc 



{\xi -X2I +ihv{Ti - T2)) 



-g/2TThv 



(5J 



It follows from (57) and (58), that in the limit of zero temperature, {h^v) ^ ^ 0, and of 
infinite size, 1/ Re -^ 0, the two-point correlation function behaves like 



r(a;i,ri;a;2,r2) 



y^pixi)p{x2] 



\XI - X2\ +ihv{Ti - T2)\^/^' 



(59) 



The latter formula is valid in the limit phv/Rc ^ 0, as well as in the limit Rc/{Phv) ^ 0. In 
(59), 9 denotes the critical exponent: 9 = 2-nhv/g, and the arguments Xi and X2, ti and r2, 
are assumed to be sufficiently close each to other. Using the notations v = ^J K/m for the 
sound velocity and p = A/g for the density of the homogeneous ideal Bose gas, we obtain 
for the critical exponent the following universal expression [11,12]: 



9 



2Txhp 



mv 



(60) 



2. The trapped Bose gas. High temperature case: /c^T ;> hv/Rc 



Let us consider the case of non-homogeneous Bose gas, which is described by the Hamiltonian 
(1) with the external potential V{x) = '^VL^x'^. In the present subsection we are mainly 
following the content of the paper [20]. Let us consider the non-homogeneous Stringari 
equation (48) and write it again, for convenience: 



-^d'^G{x,T;x',T') + d^\^\^l- —jdccG{x,T;x',T' 



h^v"^ 



5{x — x')5{ 



T — T 



(61) 
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We consider (61) for the arguments (x, r) G [—-Re, -Re] x [0,/9] with the periodic boundary 
condition only with respect to r (contrary to equation (53), 5{x — x') is a usual Dirac's 
5-function with a support at the point x' G M). The Green function satisfying (61) can be 
written as a formal Fourier series: 



G{x, r; x', r') = ^Y. ^^"^"""'^^^(a^, a:'), 



(62) 



where oj = {2n/P)l, I G Z. The spectral density G^{x,x') in (62) is then governed by the 
equation 



UJ 



h?v'' 



«"'-^') + s((^-|!)l''"<^-^'') = s4^^(^-^')- 



(63) 



Solution of equation (63) can be obtained in terms of the Legendre functions of the first 
and second kind, Pu{x/Rc) and Qy{x/ Re), which are linearly independent solutions of the 
homogeneous Legendre equation (32). As a result, we get: 



GJx^x) 



gRc 

2n?v' 



e{x-x'){Q^[^)p^ 






(64) 



where v looks as follows: 



Rr 



00" 



and e{x — x') is the sign function e(x) = sign(a;). Validity of the solution (64) can be verified 
by direct substitution of (64) into (63), where the following expression for the Wronskian of 
two linearly independent solutions P^, and Q^, [32] should be used: 



PMJyQM-QMJyPM 



[i-y 



2\-l 



and the rule of differentiation of the sign function is: (d/dx) e{x) = 26{x). 

Provided a dependence on r is neglected, an equation, which arises as a result of cal- 
culation of the correlation function accordingly to [25], looks analogously to equation (61), 
but a factor P~^ is present in its Right Hand Side instead of S{t — r'). In this case, the 
corresponding solution of the non-homogeneous equation (i. e., the fundamental solution) 
G{x; x') takes the form 

G{x\x') = -Go{x,x'), 



(65.1) 



where 



Go{x,x') 



gRc 

2h^v^ 

gRc 

{2hv) 



t[x 



In 



^''{«»(9-«»(^ 



\x-x'\ 
2Rc 



{x+x'f 
4i?2 



2Rc 



(x+x')^ 
4R2 



(65.2) 



l+x 



An explicit form for the simplest Legendre functions Po{x) = 1 and Qo{x) = | In j^^ is 
essential for obtaining the relations (65). The fundamental solution G{x]x') (65) becomes 
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equal to zero when its arguments coincide. A substitution of (65) into the representation 
(49) gives the following result for the stationary correlation function [17,18,33,34]: 



r(a;i;a;2) ~ \/ Ptf{xi)ptf{x2) exp I 



gRc 



I3{2hv) 



In 



1 + l^^^^^l - 


X^X2 

Rl 


1 \xi-X2\ 


X^X2 

R? 



(66) 



Before studying the behaviour of the correlation function, which depends on r, it should 
be noticed that solutions of the homogeneous Legendre equation (32) can be added to the 
Green function (62), since the latter respects the non-homogeneous equation. In order to 
ensure a correct asymptotical behaviour of the spectral density we are interested in at large 
|cc;|, we are free to add such a term at Ic^l 7^ and obtain the following expression: 



^iO X^"! ^ ) 



<^-^'i<iMi 



gRc 

2h^v^ 



^ , X \ / X 



Ip (^]p (^ 
2 ArJ Ar, 



(67.1) 



The Green function given by (62) and (67.1) can be represented in the form, which allows 
to study the corresponding asymptotical behaviour. In the case of strong inequality (3~^ = 
ksT ^ hv/Rc, we approximately obtain for non-zero frequencies: \uj\ ^ hv/{2R^. Let us 
take into account the following asymptotics of the Legendre functions [32,36]: 



cos( 



r(z/+l) /2\V2 1 



sm 



r(z/ + 3/2) Vvr/ sin^/^^lcos 
1 



V 



+ 



TT 



+ 0{v 



-1- 



sm 



2\'5/2 

TT/ l,,^\^ Q\^/'^ 1 COS 



[vsinOy 



+ 



TT 



(67.2) 



where e<6'<7r — e, e>0, | argz/| < 7r/2, and the notation cos6' = x/R^ is adopted. Only 
the up or the down expressions must be simultaneously chosen inside the curly brackets 
{. . . } in (67.2). Here 5 = 1 corresponds to P^, and 5 = —1 corresponds to Q^. Substituting 
(67.2) into (67.1), we determine the behaviour of G^i^x^x') at large |ci;|: 



G^{x,x') 



2hv\uj\ Vsin e sin 9' 



exp 



Re 

hv 



\u;\\9 



(68) 



When the coordinates xi, X2 are chosen to be far from the boundary of the trap, xi, 

X2 <^ Re, but at the same time the inequalities \xi — X2I <^ £i±£2 ^^^ ^^_^ — X2\ <^ Re 

are valid, the corresponding limit should be called as quasi-homogeneous. In this case, the 

function Gq{x,x') (65.2) can be approximated up to a second order (including the latter) as 

follows: 

A 



Go{x,x') 



2ffv^PTF{S) 






(69) 



Here 5" means a half-sum of the spatial arguments of the correlator, S = ^J-|^, and v is the 
sound velocity (31). In the quasi-homogeneous limit, equation (68) can be re- written in the 
following form: 

A exp {—{hv)~'^\uj\\x — x'\) 



^ uj \X, X J 



2hvpTF{S) 



\UJ\ 



(70) 
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Substitution of (69) and (70) into the series (62) leads to the answer for the Green function 
(i.e., for the correlator of the phases) (50): 

G{x,t-x',t')-—- 7^ In {2 sinh-^(|x-x'|+^^^(^-^')) )• (71) 



Therefore, the Green function (52) takes the following form at /? ^ ^ hv/Rc'. 



|smhj|^(|xi -X2I +inv{n -r2))|V^(*'' 
where the critical exponent 9{S) depends now only on the half-sum of the coordinates S: 

mv 

The result (72), which is valid for the spatially non-homogeneous case, is in a correspondence 
with the estimation (57) obtained above for the homogeneous Bose gas. Therefore, the 
estimation (72) is also concerned with validity of the appropriate condition that the size of 
the domain occupied by the Bose condensate grows faster than inverse temperature, i.e., 
with the condition hPv/Rc — > 0. 

The relation (72) can be simplified for two important limiting cases. Provided the con- 
dition 

1«^^«^ (74) 

npv npv 

is fulfilled in the quasi-homogeneous case, we obtain from (72) that the correlation function 
decays exponentially: 

T{xi,Ti;x2,T2) ~ ^/p{xi)p{x2) exp(-——\\xi -X2\ +ihv{Ti -ra)!), 

, ^^^^ (75) 



2ph^V^PTF{S)' 

The correlation length C,{S) is defined by the relation (75), which depends now on the half- 
sum of the coordinates: 

as) = ^eiS) = ^iMli^. (76) 

TT m 

In an opposite limit, 

' ' ' , ^-^ « 1 « T^, 77 

npv p npv 

the asymptotics of r(xi, ri; X2, T2) takes the following form: 



^Jp{xl)p{x2) 
\xi -X2I +ihv{Ti - T2) 



r(a:i, ri; X2, T2) ~ j^Z ^ 1 , ,fe../_ ^\\i/e{s) - ("J 



The obtained asymptotics (78) is analogous to the estimation (59), which characterizes the 
spatially homogeneous Bose gas. But the critical exponent d{S) (73) differs from 9 (60), 
since it depends on the spatial coordinates. 
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3. The trapped Bose gas. Low temperature case: ksT <C hv/Rc 

Let us pass to another case which also admits investigation of the asymptotical behaviour 
of the two-point correlator r(a;i, Ti; X2, T2). As in the previous subsection, we begin with the 
non-homogeneous Stringari equation (48), which can be written in the following form: 

d'^G{x, t; x', t') + — %/i?,) (^ (^1 - -^ j d(^^/R^)G{x, r; x' , r') j 

where the notation a = Rc/{hv) is introduced. The asymptotical behaviour can be inves- 
tigated in two cases, and these cases may be characterized in terms of a: P/a <^ 1 (the 
previous subsection) and P/a ^ 1 (see below). The functions 



n+lPr.{ 



X , 

., n>0, 



where Pn{x/Rc) are the Legendre polynomials, constitute a complete orthonormal system in 
the space L2 [—Re, Re]- This fact allows to obtain the following representation for the Green 
function G{x, r; x', r') in the form of the generalized double Fourier series: 

«(--'-'.^(^)EE^-..(^)-»(i)^--. (-) 

In (80), as well as in (54), (62), the notation ^ denotes a sum over the Bose frequencies 
ijj = {2n/P)l, I G Z, and the following notation for the energy levels (33) is adopted: 



E^ = hnf^""^'^ = ^/^^^. (81) 

Let us note that a transition from (54) to (80) takes place provided the wave functions and 
the dispertion relation are appropriately substituted. 

After summation over the frequencies and after regularization consisting in neglection of 
the term corresponding to zero values of u and n, G{x, r; x', r') (80) takes the form: 

G{x^ t; x', t') 

^ i-i 00 ^^ri / 27rAr 7\ ^ 00 -, /r-* / 

- {^9_\ \f£Y sr ^""^^^^-^ ^^\- ^ + V2 /^x fx^ 

\pRj[\27rJ ^ P ^ 2 ^ En ""KrJ^KRc 

1=1 n=l 

X ('coth('^K)cosh(E„Ar) -sinh(E„Ar)) , (82) 

where At = \t — t'\. The obtained representation (82) admits an investigation for two 
limiting cases: P/a <ti 1 (this case agrees with the estimation (75) of the previous subsection) 
and P/a ^ 1. 

Indeed, putting r = r' in the case P/a ^ 1, we obtain [35] the following relation: 

r^f I \ 9P 9 v^ n + 1/2 ^ ( x \ ( x' 

n=l " 

gRc ( \x — x'\ xx'\ gRc /I , P"^ \ 

= 2^^^H^ + ^^-^J+^^(2-^"2-^. (83) 
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Let us observe that Right Hand Side of (83) respects the non-homogeneous equation 

h'^Mi' - J)9w«..G(x.r;.'.r)) = ^«^) - i). (84,1) 

Provided the spatial arguments in the relation (83) are equated, one can obtain the following 
equality: 

G{x, r; X, r) + G{x', r; x', r) 



-|^-'[(-S)(-f)]-ii.(i---^)^ («-) 

Direct substitution demonstrates that (84.2) satisfies an equation, where Left Hand Side is 
the same as in (84.1), while only the constant term ^^ is present in Right Hand Side. 
Therefore, a subtraction of Right Hand Side of (84.2) from Right Hand Side of (83) results 
exactly in the fundamental solution G{x;x') (65). The Green function G{x;x') becomes 
equal to zero at a; = x', and it respects an equation of the type of (84.1) but with %S{x — x') 
in Right Hand Side. 

With the help of the series (83), we may bring the exponent in (49) into the form: 

(^)E^(-.(9--"(|))^- '-' 



ra=l 



As it is seen from (83) and (84.2), the relation (85) is nothing else but the fundamental 
solution G{x]x') (65) taken with an opposite sign. Thus, the representation (83), being 
substituted into (49), leads to the formula for the "equal-time" correlator r(a;i,T;a;2,T), 
which has the same form as Right Hand Side of (66). Taking into account the quasi- 
homogeneity conditions, which were discussed in the previous subsection, and recalling the 
corresponding estimation (69), we obtain for G{x]x'): 

In its turn, the latter formula leads to the estimation: 

T{xi,T;x2,r) ~ y/pTF{xi)pTF{x2) exp (-TTmki -^slj, (86) 

where C,{S) is the correlation length (76). As a result, equations (75) and (86) demonstrate an 
agreement of the estimations based on two different representations for the Green function 
G{x,t;x',t'): the first one is in the form of the trigonometric Fourier series (62) (where 
either (64) or (67) is used to express Gi_^{x,x')), and the second one is in the form of the 
series (82), which runs over the principle quantum numbers (obtained from the generalized 
double Fourier series (80)). 

Now let us turn to the case P/a ^ 1, where PEn ^ 1, Vn = 1, 2, . . . In other words, let 
us suppose that ksT <^ En and, so, ksT <^ Kl. Then, from (82) one obtains: 



G{x,t;x,t) = -^ 



1 Ar\2 1 

2 ~ ~]f) ~ T2i 
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Let us note that a difference between two neighbouring energy levels (81) can be esti- 
mated. After some appropriate series expansions, which are valid at n > 1, one obtains: 



E^, 



n+l 



E ^ — 

a 

1 

1 + -^ 



1 

1 



1 



n + l)2 
55 



+ 



7 



1 



128 {n + 1)^'" 



5n^ 



+ 



4n3 128n4 



Right Hand Side of (88) demonstrates that the levels (81) can approximately be treated 
as equi-distant provided the inverse powers of n are neglected in (88), say, for the values 
n > no = 10. In its turn, the following series expansion is valid: 



n+1/2 
A/n(n + 1) 



1 



1 



5n^ 



15 



m-^ 



128n4 



^9) 



It is remarkable that the terms ~ n~^ are absent both in (88) and (89). Let us remind that 
leading asymptotical estimations, which are obtainable with so-called logarithmic accuracy, 
are important for physical applications. As it will be clear below, the problem at hands just 
admits an estimation with leading logarithmic accuracy. In this case, the inverse powers 
of n can be omitted with the same accuracy in (88) and (89) at n > uq = 10. For such 
approximation, the energy levels (81) turn out to be equi-distant, while the corresponding 
ratio " in (87) becomes equal to unity. Convergency of the series (87) is not affected 

in this situation, while the term omitted can be estimated. 
Let us consider the exponent in (87): 

r-, rAr Ar/ 1\ Ar / 1 

Jnin +1) = (n + -) — 

^ ^ 'a a \ 2/ a \8n 



1 



5 



16n2 128n3 



(90) 



The second term in (90) can also be neglected in the exponent at n > no, provided At /a <^ 1. 
This means that the series entering into (87) can approximately be written as follows: 



no 



n + 1/2 ^/x\^/x'\ ( I--, rAr 



+ e 



-AT/{2a) 



E P- 



n=no+l 



X 
Rr 



P. 



-Ar/a 



(91) 



where no is the number, which is fixed (its specific value is forbidden to go to infinity). The 
correction, say c, omitted in the representation (91) can be estimated: 



X e 



n? 



^AT/(2a) 



a no 



E \M^y4^ 



-Ar/a 



n=no+l 

It can also be demonstrated that absolute value of the first term in (91) does not exceed 
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2no. Using (91), let us put (87) in the following form: 



G{x,T-x',r')^^ 


r/1 At\2 1i 

IV2 /3 J I2J 




2^^^i V"(" + l) V ^ ^ ' « 


xP„( 







exp(-(n+i)-^) 



00 / 

n=0 '^ '^ 

where t = exp(— Ar/a). As it can be seen from (80), the equation obtained (92) is valid for 
the case when r and r' are close either to zero or to P, as well as for the case when only r 
or t' is close to p. Besides, we assume that r 7^ r' in order to keep convergency of (92). 
Using the known series [35] 

°° 4 1 

V rP„(cos^i)P„(cos^2) = Mi^), < t < 1, (93) 

^^-^ n Ua.u_ 

n=0 + 



u+ = ^l-2tcos{^i+^2)+t\ u- = yi-2tcos(^i -^2)+t2, 



K 



(94) 



M^ + «_ ' 



where K is a complete elliptic integral of the first kind, one can estimate the approximate 
representation for the Green function (92). Indeed, let us put cos^i = x/Rc <^ 1 and 
cos ^2 = x' / Re ^ 1. Then the following estimations for m+ and U- can be obtained: 

«,^l + t-^i^±^^2-^-i^±^ (95 1) 

''+^ + * 1 + t 2RI ^ a mi ' ^^^-^^ 



//, ,n (a; — x')^\^/2 \\x — x'\+ihv/\T\ /^^ <^x 

((1 - t)2 + t^-^ j ^^ ^ ^^«., (95.2) 



where it is assumed that lS.r ja <^ 1 and 



t 1 /^ Ar 



1 + t 2 V 2a 

Provided the terms of second order smallness are neglected, an estimation /t ~ 1 — m* 
arises for k, (94), where m* implies the corresponding approximate value of m_ given by Right 
Hand Side of (95.2). When k ~ 1, several leading terms of the asymptotical expansion of 
the function 'K.{k) [35] can be written down: 

K(fi;) ^K(l-M,)^^((— + l)ln— -1), m, < 1. (96) 

The estimation presented (96) does not contain the terms ~ (m*)^, and the other terms 
which contain higher powers of m* are absent. This is due to the fact that under the imposed 
condition of quasi-homogeneity, the value m* (95.2) is treated as a quantity of first order 
smallness. Therefore, the corresponding value of the argument k, i. e., k ~ 1 — m^,, is written 
neglecting the contributions of second order smallness. 
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Let us point out that the first and the second terms in G{x, r; x', r') (92) are not small, 
because the number uq can occur to be large. Besides, the inequality 13 /a ^ 1 implies that 
gf3/Rc ^ g/{hv). However, provided a smallness of the quantities x/Rc, x'/Rc and Ar/a is 
taken into account, it is seen that the first two terms in (92) are less important in comparison 
to the third one, which can be logarithmically large at sufficiently small m*. It is why, we 
neglect the first two terms and write down the leading contribution to the Green function 
(92) with the logarithmic accuracy: 

G{x, t; x', t') 

^ (^\ fl + ^-^ ) In '-^ 

\AnhvJ \ \\x — x'\+ ihv{T — t')\/ \\x — x'\ + ihv{T — t')\ 

^ ^ In—. (97) 



^ 1 


1 1,1 


<— . 


K — < — In — 


^ M* 


M* M* M* 



27!-hv/ M* M^, 

Here it is assumed that m* is given by (95.2), and the following conditions of validity of the 
logarithmic estimation are respected: 

(98) 

Substitution of (97) into (52) gives us the following estimation for the two-point correlator 

T{xi,Ti;x2,T2): 

IFi - X2\ +tnv[Ti - T2)\'-'" 

In (99), the notation 9 for the critical exponent is introduced: 

- 27rhv 

e = M,. (100) 

The critical exponent 9 depends on m* (95.2), and therefore it is a function of differences of 
the coordinates. Apart from the inequalities (98), the following estimations can be obtained 
to characterize the relations (99), (100): 

1 i?c 13 

< ^ 7. -7—, ^ < 



hv hv\\x — x'\+ ihvij — t')\ \\x — x'\+ ihvij — t')\ 

1 /? /? 



hv Re \\x — x'\+ ihvij — t')\ 

The estimation obtained (99) (together with the critical exponent 9 (100)) constitutes 
the main result of the present subsection devoted to the case given by fc^T <^ hv/ Re. From 
a comparison with the spatially homogeneous Bose gas, one can see that the derivation of 
the estimate (99) is just analogous to a transition from the relations (56), (58) to the final 
asymptotics (59). Then, validity of the corresponding limit Rc/{hl3v) — > means that the 
result (99) is also concerned with the fact that the condensate boundary Re increases slower 
than the inverse temperature. 

Let us note that due to (98), a specific value of Uq can be related to the size of the trap 
Re- at a fixed range of deviations between the spatial coordinates x and x\ increasing of 
Re leads to an increasing of an upper bound for admissible values of no. However, due to 
(98), the estimation obtained for G{x,t;x',t') (97) does not depend explicitly on a specific 
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choice of the number no. On the other hand, at sufficiently large values of n, the following 
asymptotics for the Legendre polynomials P„ is valid [35]: 



P„(cos^) 



vrn sin {} 



cos 



n + -]'& 

2/ 4 



0(n-^/^), < ^ < 



TT. 



(101) 



Let us assume that the number no is large enough to substitute (101) into the series 



n=l 



X 



which is a part of the representation (92) (see also (91)), in order to obtain its estimation 
with the logarithmic accuracy. Eventually, the following result takes place for G{x, r; x', r') 
(92): 

G(^> -; ^\ r') - (t^) f 1 + ^) In ,., .., ,^L.._ _., , (102) 






+ ihv{T — t')\ 



where S = ^^^^. In the limit 1/Rc -^ 0, the total coefficient in front of the logarithm in 
(102) takes the value —1/9, where the critical exponent 9 is defined like in (59), (60). 

The usage of the asymptotics (101) implies that the eigen-functions are approximated 
by the base consisting of the plane waves, which correspond to an almost homogeneous 
Bose gas. Therefore, a replacement of the asymptotics (97) by the asymptotics (102) can 
be explained as a transition, at increasing of Re, to smaller scales characterized by small 
ratios x/Rc, x'/Rc and \x — x'\/Rc (the quasi-homogeneity condition). Then, the result (99) 
together with the critical exponent (100) demonstrate an effect of finitness of the size of the 
domain occupied by the spatially non-homogeneous Bose gas. This follows just from the 
employment of the Legendre polynomials as the base of one-particle states. 

§5. Conclusion 

The model considered in this paper describes a spatially non-homogeneous one-dimensional 
Bose gas with a weak repulsive coupling placed into an external harmonic potential. The 
paper deals with an application of the functional integration approach to the calculation 
of the two-point thermal correlation function of the non-homogeneous Bose gas. The tem- 
peratures that are low enough for the quasi-condensate to be created in the Bose system 
in question (see §2) are studied. The functional integral representation for the two-point 
correlation function in question is estimated by means of the stationary phase approximation 
in the way proposed in [25]. The main results are obtained for the case when the size of 
the domain occupied by the quasi-condensate increases, while the temperature of the system 
goes to zero. It is demonstrated that the behaviour of the correlation function near zero 
temperature has a power-like dependence, and it is governed by the critical exponent. In 
contrast with the case of spatial homogeneity of the Bose gas, the presence of the external 
potential is manifested in the non-homogeneity of the critical exponent. The latter turns 
out to become a function of the same spatial and thermal arguments as the correlator itself. 
The dependence of the critical exponent on these spatial arguments is in correspondence 
with the limiting behaviour of the ratio of the size of the trap to the inverse temperature 
provided both of them are increasing. 
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